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PACS 85 . 65 . +h - Molecular electronic devices 

PACS 71.70.Ej - Spin-orbit coupling, Zeeman and Stark splitting, Jahn- Teller effect 
PACS 85.85.+j - Micro- and nano-electromechanical systems (MEMS/NEMS) and devices 

Abstract. - We calculate the electronic transport through a molecular dimer, in which an excess 
electron is delocalized over equivalent monomers, which can be locally distorted. In this system 
the Born-Oppenheimer approximation breaks down resulting in quantum entanglement of the 
mechanical and electronic motion. We show that pseudo Jahn- Teller (pJT) dynamics of the 
molecule gives rise to conductance peaks that indicate this violation. Their magnitude, sign 
and position sharply depend on the electro-mechanical properties of the molecule, which can be 
varied in recently developed three-terminal junctions with mechanical control. The predicted 
effect depends crucially on the degree of intramolecular delocalization of the excess electron, a 
parameter which is also of fundamental importance in physical chemistry. 
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• Introduction. — Nano-electromechanical devices 
(NEMS) electrically detect and control mechanical mo- 
tion with great precision [1] and can be constructed in 
various nanostructurcs, including macromolecules such as 
suspended carbon nanotubes [2,3]. Nowadays even nano- 
meter sized molecules are within reach of experimental in- 
vestigation. Successful three-terminal transport measure- 
ihents [4] have been reported, detecting the quantized vi- 
brational [5,6], spin [7] and magnetic [8,9] excitations of a 
single molecule. The quantum-limited operation of NEMS 
is thus a starting point, rather than a goal in the single- 
molecule regime. More challenging is achieving control 
over such devices. Recently, electrical three-terminal de- 
vices have been demonstrated with additional mechanical 
control [10, 11]. Here the size of the nanogap in which the 
molecule is embedded can be adjusted with sub- Angstrom 
precision, thereby changing the capacitivc and resistive 
coupling of the molecule to the source and drain elec- 
trodes, as well as the intrinsic molecular properties. The 
interesting question arises how single-molecule quantum 
states involving electronic and mechanical degrees of free- 
dom may be detected and controlled in such transport ex- 
periments. This has been addressed in several theoretical 
studies, e.g. [12-15]. Fundamental to nearly all of these 
works is the adiabatic Born-Oppcnhcimer (BO) approx- 



imation, where one separates the timescales of the fast 
electronic motion from the slow dynamics of the nuclei. 
The transport is then governed by the Franck-Condon 
(FC) principle, where the tunneling of an electron onto 
the molecule changes the electron number N N + 1, 
which induces a transition between electronic states e e' 
and the initial vibrational state Xe is projected onto a vi- 
brational state x'e' of the molecule in the final electronic 
state. The amplitude for this process factorizcs and is 
proportional to the overlap integrals of the mechanical 
wave functions {Xe\Xe') which in general is non-zero and 
strongly depends on both of the vibrational states. There 
are thus no selection rules in contrast to spin-related tun- 
neling. By proper choice and design of the vibrational 
properties of molecular transistors (number of modes, adi- 
abatic potential landscapes, etc.), the FC-effect may thus 
be exploited to induce non-trivial vibrational states and 
interesting non-linear transport characteristics relevant for 
possible electro-mechanical sensing applications. 
A novel aspect of molecule-based NEMS is they may dis- 
play strong vibronic effects (distinct from vibrational) due 
to the non-trivial coupled quantum dynamics of the elec- 
tronic and nuclear degrees of freedom, for a review see [16]. 
Here the Born-Oppenheimer separation of the time-scales 
for the nuclear and electron motion breaks down. The 
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system is only adequately described by so-called vibronic 
states, in which the quantum entanglement renders the 
concept of electronic and nuclear motion meaningless. The 
most prominent and well studied vibronic effect is the dy- 
namical Jahn- Teller (JT) effect, which occurs in molecules 
where the electronic ground state is degenerate due to 
a high spatial symmetry of the static nuclear framework 
of the (non-linear) molecule. In such systems, there al- 
ways exists [17] a vibrational coordinate along which a 
static distortion will lower the molecular symmetry and 
lift the degeneracy. However, in a single molecule, the 
distortions are dynamical and the electronic degeneracy 
is transformed into a vibronic degeneracy i.e. of the 
quantum-mechanical molecular eigenstates. Recently, the 
selection rules encoded in these molecular eigenstates (re- 
lated to the high symmetry) were predicted to block elec- 
tron transport through a JT active molecule [18]. An 
important question now is how to distinguish such vi- 
bronic blockade from spin- [19], magnetic [20] or Franck- 
Condon [14,21,22] blockade effects. Another issue is that 
the BO-approximation breaks down even when electronic 
levels only come close in energy (on the order of a few 
vibrational energy spacings). This is referred to as the 
dynamical pscudo Jahn-Teller (pJT) effect and occurs in 
many molecular systems [16]. A generic problem where 
it occurs is in determining the degree of delocalization 
of an excess electron in molecular mixed-valence com- 
pounds [23] , which is fundamental to the classification by 
the Robin-Day (RD) scheme in physical chemistry [24]. 
The dynamical pJT effect is relevant in any system where 
electrons become delocalized over multiple similar centers, 
while local vibrational modes couple to the charge on the 
centers. Molecular dimers, being of recent experimental 
interest [6,25], constitute a basic system which may ex- 
hibit this effect. Such a system may be called "molecular 
double quantum dot" . However, due to the vibrations and 
the pJT effect its transport properties dramatically devi- 
ate from semiconductor double quantum dots. More gen- 
erally, such mixed-valence effects are important in func- 
tional nanosystems such as metal-organic supramolecular 
arrays [26] and polyoxometalates, see e.g. [27]. 
In this letter, we predict transport signatures of the pJT 
dynamics, of a single molecular dimer transistor, which 
markedly differ from those due to the Franck-Condon ef- 
fect, allowing the breakdown of the Born-Oppenheimer 
principle to be identified experimentally. We find charac- 
teristic non-linear conductance peaks with a sharp depen- 
dence of their position, magnitude and sign on the electro- 
mechanical parameters of the molecule. 

Model. — We consider a dimer molecule consisting 
of two identical monomers, labelled hy i = 1,2. Each 
monomer can vibrate along its individual totally symmet- 
ric ("breathing") mode Qi about the potential minimum 
at Qi = with frequency lo. Each monomer also acco- 
modates one electronic orbital state, \ia), for an excess 
electron with spin projection a that can tunnel to the 



other monomer with amplitude t via a mechanically stiff 
bridging ligand. It thereby significantly distorts the occu- 




Fig. 1: (Color online) (a) Sketch: dimer molecule trapped in a 
nanogap between two voltage biased electrodes, Vt — — /J-r, 
and capacitively coupled to a back-gate (not shown) at volt- 
age Vg, which shifts the effective molecular energy levels. The 
monomers (blue) can vibrate along the local totally symmet- 
ric breathing mode and are connected by a mechanically stiff 
bridge (not shown), (b-c) Adiabatic potentials Wg and We 
for the symmetry-breaking molecular distortion Q_ (full lines) 
and harmonic expansions around Q- = (dashed lines) for (b) 
weak (A = 0.7, t — 2.15a;) and (c) intermediate strength of the 
pJT effect (A = 1.93, t = 2.15a;). 

pied monomer along coordinate Qi due to a change of the 
bond- lengths (c.f. Fig. 1(a)). The resulting shift of the 
potential minimum is \/2X, expressed in units of the zero- 
point motion energy of the vibration of the undistorted 
monomers. Thus A is the dimensionless electron-vibration 
coupling. This model applies, for example, to a mixed- 
valence molecule, where the monomers are metal-ions 
(zero spin) with a surrounding shell of ligand atoms. Be- 
cause of its coupling to the presence of the hopping excess 
electron, the local distortion of the monomer is dragged 
along by the electron and becomes similarly delocalized 
over the dimer. This results in coherent electro-mechanical 
motion and the breakdown of the Born-Oppenheimer sep- 
aration. The central quantity controlling the character of 
the molecular states is the delocalization energy t relative 
to the coupling Xoj to the localized distortion. We note 
that in ref. [28], the opposite case of stiff monomers and 
a distorted bridge was considered, requiring only a single- 
mode to be considered, whereas here we account for two 
vibrational modes and their interplay. Also, below we con- 
sider transport up to high bias Vb ^ At, in contrast to [29]. 
When the size of the nanogap is varied, which is possible 
in a mechanically controllable break junction, the intra- 
molecular hopping, t, will change. This allows for an in 
situ change of the character of the molecular state. In 
principle, changes in cu and A may also occur, but they 
do not alter the results qualitatively and can safely be ne- 
glected. 

We assume charging effects (Coulomb blockade) to be 
strong enough that only two molecular charge states par- 
ticipate in transport processes, which we label by the 
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number of excess electrons on the molecule N = 0,1. 
The Hamiltonian for the molecule in charge state N, 
written in the molecular vibrational coordinates Q± = 
{Q2 ± Qi)/V2, then reads 

= H'^ -XLuQ+ + XLuQ^{ni~h2) 

+tY,idLd2a + h.c.) (2) 

where d]^ creates an electron in state |icr) and fii = 
Xio- d\„dic is the occupation operator. The symmetric 
coordinate, Q+j corresponds to the monomer vibrating 
in phase i.e. the molecular breathing mode where the 
molecule as a whole changes its size. It couples to the 
total excess charge, oc A^, of the molecule, resulting in 
a simple shift of the potential surface along by an 
amount A (linear term in eq. (2)). The Franck-Condon 
(FC) transport effects resulting from this type of cou- 
pling have been calculated by several groups [13-15, 21] 
and found experimentally [5,6,30]. In contrast, the anti- 
symmetric mode, Q-, corresponds to the monomers vi- 
brating with opposite phase. If an excess electron is 
present, this molecular shape distortion couples to the 
internal charge imbalance hi — 77.2 • Due to the intra- 
molecular tunneling, t, the Hamiltonian (2) mixes elec- 
tronic and vibrational states of the mode Q- prohibiting 
a factorization of the molecular wave function into a. 
vibrational and an electronic part. We thus need vibronic 
states, |m_,o-) = |Xm_)|l''') + IXm_)|2o'), to describe the 
excess electron and the pJT-active mode. Here m_ de- 
notes the vibronic quantum number for the joint electron- 
vibration (Q_) system: distinguishing these systems is 
fundamentally impossible due to the quantum coherence. 
Finding the vibrational coefficients \Xm_ ) in the molecular 
vibronic eigenstates requires diagonalization of the Hamil- 
tonian (2), which has to be done numerically. Despite the 
breakdown of the BO-approximation, it is instructive, e.g. 
for interpreting the energy spectrum, to consider the adi- 
abatic potentials for the Q- vibrations, obtained by ne- 
glecting the nuclear kinetic energy operator (P_ 0) in 
eq. (2) and to find the electronic eigenstates as function of 
while neglecting Q+. The resulting electronic energies 
are the ground (g) and excited (e) adiabatic Q_-potential 
for the vibrations, Wg^e{Q-) = ^(^Q- T \/{^^Q-Y + 
which are sketched in fig. l(b-c). However, in the sec- 
tion Results we will discuss when the breakdown of the 
BO-approximation renders these potentials meaningless. 

Transport. — The transport setup, sketched in fig. 1 
(a), is described by the Hamiltonian H = Hm + Hrcs + Ht- 
Here the molecular part Hm is specified by — aVgN 
for N = 0, 1 excess electrons where the gate voltage Vg 
linearly shifts the addition energies (proportional to the 
gate coupling a). The electrodes r ^ L,R arc described 
by iJros = Y.r,k,a i'^k - l-h) clf,^Crka and are kept at tem- 



perature T and electrochemical potentials = M i Vf,/2 
i.e. Vb is the bias voltage. They are coupled to the 
molecule by Ht = J2r,i,k,cr'^rdlaCrka + h.c, where cl^^ 
creates an electron of spin a in state k in electrode r and 
is the amplitude for tunneling between the electrode r 
and monomer i — 1,2. We consider a linear arrangement 
as sketched in fig. 1(a) by assuming T| = = and, 
for simplicity, symmetric coupling Tl = = yjT / {^irp). 
Here r denotes the tunneling rate and p is the electrode 
density of states. In many experiments [5-9, 30] the 
transport is dominated by single-electron tunneling. We 
therefore focus on the corresponding range of applied 
voltages and temperature. We calculate the non- 
equilibrium stationary state occupations of the molecular 
states and the transport current both up to second 
order in iJx using a standard kinetic (master) equation, 
see e.g. [15]. We have checked by explicit calculation 
that the off-diagonal elements of the density matrix 
are not crucial. Rather it is more important to take a 
large number of states into account (order of 1000 states). 

Results. — In order to appreciate the breakdown of 
the BO separation, we first discuss a case where it has 
approximate validity i.e. the effect of the distortion of the 
monomers is sufficiently weak. In fig. 2(a) we show the 
differential conductance as function of the applied volt- 
ages. Many excitations appear, involving the and/or 
the mode, which are separated in bias voltages by mul- 
tiples of 2u} (due to symmetric biasing). The first of these 
excitations starts out at the marker (i) in fig. 2(a). How- 
ever, in contrast to usual FC transport spectra [13-15], 
the Q- excitations are weakened within the noticeable 
gap of At (due to exponentially suppressed overlap inte- 
grals of classically forbidden transitions), and enhanced 
conductance peaks delimit the upper boundary of the gap 
(starting out from marker (ii)). One can thus directly esti- 
mate the strength of the delocalization of the excess elec- 
tron. Furthermore, the excitations spaced in VJ, roughly 
by 2u} also have a detailed substructure of a dense series 
of conductance peaks, for instance, along the right edge 
of the transport region, terminating at the marker (iii). 
The latter correspond to tunnel processes between Q- ex- 
citations of the = and = 1 ground potential with 
the same vibrational quantum number m_ = 0, 1, 2, . . .. 
For A^w ^ t, the potential Wg{Q-) is approximately har- 
monic, but with reduced frequency LOgjuj ~ 1 — ujX"^ /t due 
to the pJT interaction, see fig. 1(b). Therefore the res- 
onances corresponding to different to_ occur at slightly 
different positions [15]. The equidistant energy spacings 
correspond to uj~u)g ~ {Xuj)'^ /t <C uj. Similarly, above the 
gap near marker (ii), a dense series of conductance peaks 
with negative Vg-dependence indicates that the upper adi- 
abatic potential has a higher frequency We/cj = 1 +10X^/1. 
Using the gap and features at (i)-(iii) in fig. 2(a) one can 
thus estimate A, t and oj from the transport data. For 
larger values of A, the adiabatic potentials additionally 
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Fig. 2: (Color online) dl/dVt (F = 2.5 • IQ-^c^, T = 4 ■ IQ-^o;) 
for (a) weak pJT mixing (A = 0.7, t = 2.15[i)) (b) moder- 
ate pJT mixing (A = 1.93, t — 2.15cj). Inset: high contrast, 
dashed black line marking the Vb — Vg trace taken in fig. 3(a). 
For convenience the gate voltage is defined such that Vg = 
corresponds to the charge degeneracy point. Due to symmet- 
ric biasing, energy scales appear at twice the separation on the 
voltage axis. 




become anliarmonic due to the pJT interaction resulting 
in markedly non-equidistant spacing of the dense series of 
conductance peaks. We note that, interestingly, in mixed- 
valence molecules with magnetic ions this vibrational an- 
harmonicity correlates with the relative orientations of the 
ionic spins [31]. 

The most dramatic effect of the pJT interaction is the 
breakdown of the BO separation for stronger coupling A. 
This occurs when excited states of the two adiabatic po- 
tentials We and Wg come close in energy. These anti- 
crossings can be observed in the transport at high bias 
voltages, Vb > 4i (c.f. fig. 1(c)), where the states of the 
excited adiabatic potential become accessible as well. At a 
first glance, the transport spectrum in this case, shown in 



Fig. 3: (Color online) Signature of the pJT effect: anti- 
crossings as the intramolecular delocalization t is varied due 
to a mechanical change of the nanogap size, (a) Evolution 
with t of the dl jdVb trace along the line in the Vb, Vg plane 
marked in the inset of fig. 2(b). (b) Evolution with t of ener- 
gies of the vibronic states for A'' = 1 (A = 1.93), with the har- 
monic (5+-vibration energies subtracted. The green (dashed) 
/ red (solid) color indicates positive / negative parity. The 
anti-crossing in the transport in (a) corresponds to the marked 
anti-crossing around t ~ 2.3lj and El_ ~ buj. (c) Evolution 
with t of the vibrational parts xln_{Q-) for = 9 and 10, 
respectively. Due to the symmetry of the dimer, the vibronic 
state has definite molecular parity tt = ±, reflected by the 
property Xm_(Q-) = ■^xli.A-Q-)- 
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fig. 2(b), seems inextricably complex. However, clear sig- 
natures of the pJT effect are revealed when the nano-gap 
size is varied and the intramolecular hopping t changes 
while the mechanical properties of the monomers A, to re- 
main fixed. In fig. 3(a) wc show the evolution of the dl/dVb 
trace taken along the dashed black line in the inset of 
fig. 2(b), as t is varied. Experimentally, such data can be 
collected with techniques described in [10,11]. Among the 
most pronounced dl/dVb resonances, those with a weak t- 
dcpendcncc one would assign to highly excited vibrations 
in the lower adiabatic electronic state, Wg, whereas those 
with a strong linear t-dependence would correspond to the 
lowest excitations in the upper adiabatic electronic state. 
We- The main difference in t-dcpendence stems from the 
gap, 2t, separating the adiabatic potentials (c.f. 1(c)). 
This distinction is completely lost at the anti-crossings vis- 
ible in fig. 3(a). The conductance anti-crossing at t w 2.3a; 
maps out the corresponding anti-crossing in the evolution 
of the vibronic energy spectrum with t which is shown in 
fig. 3(b). Strikingly, the sign of the conductance of the two 
anti-crossing transport resonances is different. This di- 
rectly relates to the large difference in the kinetic energy 
of the nuclear motion of the two anti-crossing adiabatic 
electronic states. In fig. 3(c) we show the real space repre- 
sentation of vibrational components of the vibronic wave 
functions of the involved states. The weakly i-dcpendcnt 
excitation has a rapidly varying wave function leading to a 
small overlap with the vibrational groundstate of the un- 
charged molecule. The occupation of this state on average 
reduces the contributions to the current of other states, 
therefore leading to negative differential conductance [15]. 
In contrast, the strongly t-dependent excitation is more 
similar in shape to the vibrational groundstate leading to 
a much larger overlap integral and therefore positive dif- 
ferential conductance. At the anti-crossing the strong pJT 
mixing causes the components of both vibronic wave func- 
tions to rapidly vary. As a result the conductance peaks 
disappear in a narrow range of t values in the anti-crossing 
region. Note that all other resonances, which follow from 
the BO approximation and the FC-principle, smoothly 
depend on t, making the pJT effect clearly stand out. 
Strikingly, the transport anti-crossings seen in fig. 3(a) 
are replicas of one and the same anti-crossing marked in 
fig. 3(b), due to the simultaneous excitation of the Q+ 
mode. This can be seen from both the voltage distance of 
the anti-crossings and the identical <-dependence. Thus, 
interestingly, the pJT-inactive mode proliferates the vi- 
olation of the adiabatic BO separation in the transport. 
The many other anti-crossings in fig. 3(b) result in pJT 
resonances at different voltages (not shown), and compar- 
ison with these calculated levels allows one to estimate the 
parameters. More generally, the effects exemplified above 
may be expected whenever the pJT mixing is important, 
that is, for minimal separation of the adiabatic potentials 
on the order of the vibrational quanta, t ^ lu, and moder- 
ate to strong coupling to the distortion, A > 1. 



Conclusion. — In this letter we have shown for a 
minimal model exhibiting the dynamical pJT effect, that 
the breakdown of the Born-Oppenheimer separation of the 
electronic and vibrational motion in a molecular transistor 
leads to novel transport resonances which can be distin- 
guished from standard Franck-Condon effects. The com- 
bination of electrostatic gating and mechanical control is 
crucial to unravel such complex molecular transport pro- 
cesses and demonstrate electro-mechanical quantum en- 
tanglement in molecule-based NEMS. Interesting candi- 
date devices are mixed-valence molecules with a moderate 
degree of intramolecular delocalization of the excess elec- 
tron, so-called Robin-Day Class II systems. Their electron 
transport properties may shed new light on the fundamen- 
tal issue in physical chemistry of their classification by 
intramolecular charge transfer. 
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